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CN ' Abstract 

^: 

^P ■ We measure the sphaleron rate for hot, classical Yang-Mills theory on the lattice, 

in order to study its dependence on lattice spacing. By using a topological definition 
of Chern-Simons number and going to extremely fine lattices (up to /3 = 32, or lattice 
p ■ spacing a = l/{8g'^T)) we demonstrate nontrivial scaling. The topological susceptibil- 

JL . ity, converted to physical units, falls with lattice spacing on fine lattices in a way which 

Q I is consistent with linear dependence on a (the Arnold-Son- Yaffe scaling relation) and 

>-^ ■ strongly disfavors a nonzero continuum limit. We also explain some unusual behavior 

of the rate in small volumes, reported by Ambj0rn and Krasnitz. 
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1 Introduction 

Baryon number is not a conserved quantity in the standard model. Rather, because of the 
anomaly, its violation is related to the electromagnetic field strength of the SU(2) weak group 

d.J'B = Nc^^e.^c.^TrF^-'F"^ = Ng£-,E^B^ , (1) 

where Nq = 3 is the number of generations.^] The right-hand side of this equation is not 
surprisingly a total derivative, with an associated charge called the Chern-Simons number. 



J^ = J^j ^'^Jb = Ncs ^ (integer) + J^ / d'xe,,, [f^^AI - ^-U,A1A]AI) , (2) 

iVcs(ti)-iVcs(t2) = r dt f d^x^EtBt . (3) 

Chern-Simons number A'^cs has topological meaning; its change in a vacuum to vacuum 
process is the (integer) second Chern class of the gauge connection. Note that, according to 
the first equation, the total baryon number Nb need only be an integer in vacuum, when 
A'^cs is an integer. The baryon number in vacuum fixes the constant of integration in the 
definition of A'^cs- 

In vacuum, baryon number can be violated by a vacuum fiuctuation large enough to 
have a nonzero integer Chern-Simons number. The efficiency of baryon number violation 
by this mechanism is totally negligible |]l|; but at a sufficiently high temperature baryon 
number violation can proceed by thermal excitations which change A^cs, and the rate for 
such a process is not necessarily very small 0. This can have very interesting cosmological 
significance, since it complicates GUT baryogenesis mechanisms and opens the possibility of 
baryogenesis from electroweak physics alone. This motivates a more careful investigation of 
baryon number violation in the standard model at high temperatures. 

The baryon number violation rate relevant in cosmological settings can be related by a 
fiuctuation dissipation relation p, ^, |[ to the diffusion constant for Chern-Simons number, 

r.Umlim«-^^-<'>-^"'°>'\ (4) 

y^oot^oo Vt 

where the angular brackets () represent a trace over the thermal density matrix. This quan- 
tity is called the sphaleron rate for historical reasons. There has been some controversy not 
only to its size in the symmetric electroweak phase (which closely resembles pure Yang-Mills 
theory, an approximation we will make from here on) but even to its parametric depen- 
dence. On purely dimensional grounds, at high temperature it must scale as T'^, but the 
dependence on the coupling constant has been more controversial. Since the natural non- 
perturbative length scale of hot Yang-Mills theory at weak coupling is ~ l/{g^T) it was long 
believed that T oc a'^T'^. In this case, T could equal its value in classical Yang-Mills theory, 
as originally suggested by Grigoriev and Rubakov 0. 

^There is also a contribution from the hypercharge fields, but it will not be relevant here because the 
topological structure of the abelian vacuum does not permit a permanent baryon number change. 



More recently Arnold, Son, and Yaffe (ASY) have argued that while the natural 
length scale for hot, weakly coupled Yang-Mills theory is ~ l/{g'^T), the natural time scale 
is different; because of interactions between the nonperturbative infrared excitations and 
essentially perturbative but very numerous UV excitations, the time evolution of IR Yang- 
Mills fields should be overdamped and the natural time scale for their evolution should be 
~ 1/((7^T), or, restoring h, ~ l/{hg^T). The appearance of h in this expression precludes 
any simple correspondence between the classical theory and the quantum theory. They then 
argue that, since only the nonperturbative IR fields can contribute to the diffusion of Nqs, 
the correct parametric behavior for F is F oc a^T^. The analogous expression in the classical 
theory has some UV regulator serving the role of h in the quantum theory; for instance, for 
Yang-Mills theory on a lattice under the standard lattice action, the UV regulator scale is 
the inverse of the lattice spacing a, so the natural time scale should be of form l/((yf'^aT^), 
leading to F oc aa^T^ |^. Subsequently, Bodeker showed that the coefficient of the a^ law 
should contain a further logarithmic dependence on ^f^ 0, or, on the lattice, on g'^aT. 

The arguments of ASY are still considered somewhat controversial. While numerical 
simulations of classical Yang-Mills theory, supplemented with added degrees of freedom in- 
tended to serve the role of the "hard" quantum UV degrees of freedom, clearly support their 



arguments |]T0|, results for F in pure classical Yang-Mills theory on the lattice [jTT], |T^, |T3[ 
have never convincingly displayed linear scaling in lattice spacing. Furthermore, the more 
recent work of Ambj0rn and Krasnitz finds two other results which are problematic to inter- 
pret if Arnold, Son, and Yaffe are correct; they find overly strong lattice spacing dependence 
for F in small volumes, and they find unexpectedly rapid falloff for unequal time. Coulomb 



gauge fixed correlators [0 



We will not address the question of unequal time. Coulomb gauge correlators here, except 
to note that we believe such correlators should show a strong volume dependence even at 
fixed k. Since the physical volume was varying along with the lattice spacing in Ambj0rn 
and Krasnitz' results, it is difficult to disentangle these two dependencies. We leave settling 
this problem to future work, but it is our general belief that such correlators will not prove 
very useful probes of infrared dynamics. 

This work is intended to answer the other two questions about the applicability of the 
ASY picture to classical Yang-Mills theory. First, we present results for F on a much wider 
range of lattice spacings, including much smaller spacings a, than has previously been done. 
Although the results at fairly large a show weak lattice spacing dependence, as a becomes 
smaller a strong a dependence sets in, which is incompatible with a'^ scaling, but is in 
accordance with expectations if Arnold, Son, and Yaffe are correct. We also present a re- 
analysis of the behavior of F in a small, fixed volume. When care is taken to ensure that the 
physical volume really remains fixed as the lattice spacing is varied, we find F to depend on 
a slightly more weakly than in large volumes. This is also expected in the ASY picture. 

2 Expected scaling behaviors 

To explain the different proposed scaling behaviors for F, we will briefly review the ther- 
modynamics of classical Yang-Mills theory. The partition function for classical Yang-Mills 
theory is equivalent to the path integral for three dimensional, quantum Yang-Mills the- 



ory with an added adjoint scalar, as originally shown by Anibj0rn and Krasnitz |]rT|. The 
classical Yang-Mills partition function, absorbing g into the definition of the connection so 
covariant derivatives are Di = di + iAi, is 

Z = JvA'}VE^S{{D,E,r)exp{-H/T), (5) 

H = Jd'xI^-^^Ft^Ft^ + lE^E^y (6) 

The delta function enforces Gauss' law. It is convenient to rewrite it by introducing a 
Lagrange multiplier, 

S{{D,E,r) = JvA-,exp UJaI^D^eA , (7) 

where Aq has the same normalization as Ai. The integral over E is then Gaussian. Performing 
it yields the partition function 

Z = JvA^VAoexpi-H/T), (8) 

H/T = Jd'x(^-^Ft^F^^ + ^{DM''iDAory (9) 

Physically the Aq field corresponds to the time component of the connection, but if we choose 
to interpret the result as a path integral for a 3-D quantum field theory then the Aq field 
just corresponds to some massless adjoint scalar. 

The form for the partition function coincides, except for the absence of a mass term for 
the Aq field, with the path integral for the full quantum theory in the dimensional reduction 



approximation [|1J, 0, |l6l, and so it quite accurately reproduces the thermodynamics of 
infrared Yang-Mills fields. This motivates the belief that the dynamics of the full Yang- Mills 
theory will also coincide with the dynamics of the classical theory. Looking carefully at Eq. 
(P), we observe that g^ and T only enter in the combination g'^T. Now P has engineering 
dimension 4, so since T is the only dimensionful quantity in hot Yang-Mills theory at weak 
coupling, we must have P oc T'^. The "naive" scaling argument for P corresponds to requiring 
that each T carry a factor of g"^, as is motivated by the form of the partition function. 

The problem with this argument, and with the argument that the dynamics of the classi- 
cal and quantum theories should correspond, is that it neglects the effects of UV divergences. 
For thermodynamic quantities the UV divergences of the classical theory arise from a finite 
number of graphs and can be absorbed by counterterms. However, from the way that the 
Aq field arose above, the classical theory is obliged to have a zero bare mass squared for this 
field. Hence the physical mass squared will approximately equal the linear one loop divergent 
contribution, mj^ oc g^TA, with A some UV cutoff. In particular, on the lattice and using 
the conventional (Wilson) action, 1/a serves the role of A, and [p!4| p| 

ml = — — , S = 3.175911536... . (10) 

vra 



^The result in the reference differs by a factor of 5/4 because it includes a contribution from a Higgs field, 
which is absent here since we treat pure Yang-Mills theory. 



The thermodynamics of the gauge fields Ai do not care about this divergence, since their 
thermodynamics have a finite hmit as mj^ — > oo. Therefore the argument that the only 
natural length scale for the Ai fields is 1/g'^T, remains valid in the presence of UV divergences. 
However, since the dynamics of the Ai and Aq fields are intertwined, it is not at all clear 
that unequal time phenomena involving the Ai fields should be UV regulation insensitive. 

Arnold, Son, and Yaffe examined the propagator of the full quantum theory at soft 
momentum, including the UV infiuences at leading order in the coupling by including the 



hard thermal loop (HTL) self-energy contribution ||T^ . For spatial momenta p in the para- 
metric regime p -C gT they conclude that the evolution of the transverse (magnetic) degrees 
of freedom responsible for baryon number violation is overdamped, of form 

^^-^^(.).(no.e). (U) 

The perturbative treatment which gives this result holds provided that p ^ g'^T, but it 
breaks down when p is of order g'^T, the scale which we are interested in. Ignoring this 
difficulty and applying it for p ~ g'^T, we find that the time scale for significant change in A 
is ~ mjj/g^T^, which in the continuum theory is ~ l/{g'^T) and on the lattice is ~ l/{g'^aT'^). 

Of course it is problematic to apply Eq. ( pT]) beyond its range of applicability. However, 
for g sufficiently small, we may apply it at the scale p ~ g^~'^T to find that the natural time 
scale here is (1/t) ~ g^~^'^T. Since the time scale for field evolution at p ~ g^T cannot be 
faster than that at a larger value of p, Eq. (^) does place a bound on the natural time scale 
for fields with p ~ g'^T; the time scale cannot differ from the estimate ~ l/{g^T) by a nonzero 
power of g. Thus, Eq. (|rT]) is enough to ensure the ASY result up to corrections weaker 
than any power of g. A much more careful study of the HTL effective theory by Bodeker 

IB[ finds that log(l/(7) corrections do occur; however, they prove to be numerically small 



Note however that the ASY argument is a parametric treatment which relies on a large 
separation between the gT and g^T scales. Numerical results, for instance for the subleading 



contributions to the Debye screening length [^, suggest that g (or, on the lattice, g aT) 



may need to be fairly small before such a separation of scales really exists. Therefore, we 
might expect that the ASY scaling behavior only sets in at reasonably small lattice spacings 
a. This motivates the study of the sphaleron rate on very fine lattices, which we take up in 
the next section. 

3 Results: large lattices, fine spacings 

To address the question: "How does F depend on lattice spacing a when a is small?" we 
track A^cs on the lattice. We use the conventional Kogut-Susskind lattice action |2I|], which 
is the Minkowski time analog of the Wilson action for Euclidean lattice Yang-Mills theory. 
We use the same discrete implementation for the time evolution as in Ambj0rn et. al. |^ 
and almost all subsequent work. Besides the size of the lattice there is one variable, /5l, 
which is a reciprocal temperature in lattice units. At tree level it is /3l = 4:/{g'^aT), but this 
relation receives radiative corrections because the UV lattice and continuum fields behave 
differently. These are treated, for the case rrijj <ti l/a^, in p3[ and extended to larger ttid 



in Appendix B of [|10|- We use the expression from there for the one loop improved relation 
between /3l and lattice spacing a, 
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(/3l - O.Q)g'T 



(12) 



where ^ = 0.152859 . . . and the functions ^{x) and S(x) are defined in |T^; the approximation 
that the correction term is 0.6 is good to 10% for all a considered here but we will use the 
full expression. In this section the difference between the naive and improved match will 
not be important, but in the next section it is essential. Henceforth, when we refer to /3l we 
will mean the "unimproved" quantity appearing in Eq. ([T^); while the quantity f3 will mean 
/3 = 4:/{g'^aT) using the improved relation for a. 

The measurement of Chern-Simons number deserves some comment. Early work on 
Chern-Simons number diffusion [^, |Il|, ^ ^ used "naive" definitions of A'^cs in which 
the right hand side of Eq. (^ is implemented as the integral of a local operator over the 
unmodified lattice fields. This approach is not topological because such a local operator on 
the lattice is never a total derivative, as it should be for Nqs to have topological meaning. 
Because of this, such a definition of A'^cs shows diffusive behavior even when there is no true 
diffusion of baryon number occurring [^, |^ . The response to true topology change can also 
get renormalized by UV fluctuations. The latter problem gets less severe as a is reduced, 
but we find that the "spurious" diffusion per physical 4-volume grows as a~^ and is therefore 
disastrous. Therefore we should use a topological definition of A'cs- 

Technically topology is not well defined for lattice fields, but it can be well defined on a 
restricted class of lattice fields which are sufficiently "smooth" [^, ^, 



In our context 

topology is well defined for suitably small lattice spacing; in practice there is no problem 
if /3 > 6 (a < 2/3g'^T), which will be the case for almost all lattices we consider. For the 
current application, two topological means have been developed; the "slave field" method 
T^ , similar to the method of Woit p8|; and "calibrated cooling" [^, an improvement 
on the field smearing proposal of Ambj0rn and Krasnitz |T^. The "slave field" method is 
numerically efficient but noisy, which means in practice that a longer numerical evolution is 
needed to get good statistics. Therefore we will use "calibrated cooling". The philosophy 
of the method is that the topological content of the connection cannot be modified by 
small local changes; therefore we may "smear" the connection, removing UV noise which 
is responsible for the poor performance of the "naive" measurement method. After the 
smearing we measure A"cs by integrating an Ola"^) improved local operator for E^Bf, and 
we cure the slight residual error in the algorithm by "calibrating" it with occasional coolings 



all the way to vacuum; the full details can be found in [30 



The problem with the above approach is its numerical cost; measuring A"cs as a function 
of time is much more expensive than generating the Hamiltonian trajectory. This is only 
really a problem on very fine lattices, as numerical cost rises as a~^. However, in this case 
the topology changing configurations we are after are many, many lattice spacings across and 
are highly over-resolved; no topological information is lost by "blocking" the lattice and then 
measuring A'cs- That is, we construct a "blocked" lattice with B times the original lattice 
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Figure 1: Right: A^cs measured with and without 2- fold blocking on a 32^ lattice at /3 = 12 
(a = l/Sfyf^T). The smearing depth is the same, 7.5a^ = 0.83/ g^T'^ for the two curves, 
with blocking (upper curve) and without (lower curve, shifted down for clarity). Left: the 
difference between the curves, which is small and spectrally white; note the difference in 
scale between the plots. 



spacing by making the connection between neighboring points on the blocked lattice equal to 
the product of links between the same points on the unblocked lattice; the details are in [Q. 
We do this first, before applying any smearing. This introduces some white noise into Nqs 
but does not change the diffusion at all, which means that results for F will be unchanged. 
Figure |l| shows a test of blocking in which we track A'^cs for the same Hamiltonian trajectory, 
with and without blocking. The difference is small and purely spectrally white. On the other 
hand, the difference in numerical effort is enormous; without blocking, measuring A'^cs took 
5 times as much CPU time as updating the fields, while with blocking it took around 10% 
as much. We will block for all data with a < l/3g'^T in this work, which means that the 
CPU time taken to measure A'^cs is negligible for all the most numerically intensive cases. 

We should also be sure to use a large enough volume to eliminate finite volume systematics 
(effectively, to perform the ^ ^ cxd limit in Eq. ^). To do so we measure F as a function 
of L = V^^^ at a fixed lattice spacing, a = l/2g'^T. The result is shown in Figure ^ The 
dependence on L is in general accordance with the results of |T^, except at small volumes 
where ours go to zero and theirs do not because of the spurious UV contributions in their 
definition of A^cs- We will use L = 10/g'^T {N = 2.5/9) for all large volume results. 
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Figure 2: Volume dependence of the sphaleron rate in cubic volumes L on a side with 
periodic boundary conditions and a = l/2g'^T {[3 = 8). Large volume behavior is obtained 
by L = 8/g'^T. Around L = 5/g'^T the sphaleron rate falls off abruptly and it is virtually 
zero already at L = 3/g'^T. 
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3.0 


0.0023 ±0.0016 


7 


3.5 


0.029 ±0.008 
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4.0 


0.112 ±0.012 
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4.5 


0.276 ±0.015 


10 


5.0 


0.545 ± 0.035 


12 


6.0 


1.18 ±0.04 


14 


7.0 


1.40 ±0.05 


16 


8.0 


1.61 ±0.07 


20 


10.0 


1.68 ±0.03 



Table 1: Volume dependence of the sphaleron rate for an a = l/2g'^T lattice. The volume 
dependence is very strong at L = 4:/g'^T but is consistent with zero above L = 8/g'^T. 




Figure 3: Chern-Simons number diffusion constant plotted against lattice spacing a, and a 
linear fit for the last 5 points. Around a ~ 0.5/g'^T the dependence is somewhat weak, but 
at larger a a rapid falloff is evident. 
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0.94 ±0.03 
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202000 


0.626 ±0.030 


0.67 ±0.03 



Table 2: Sphaleron rate F as a function of lattice spacing, f3 = 4/g'^aT. On the finer lattices 
an approximately linear dependence on a (or 1//9) becomes apparent. For those who prefer 
it we also include F converting to physical units using the naive relation between /9l and a. 
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Figure 4: T /a against a. The straight hne fit assumes the ASY scahng argument, the curved 
fit includes Bodeker's logarithmic correction to this scaling argument; each are based on the 
5 highest points. 

Our results at large volumes are presented in Table ^ and plotted in Figure ^. These 
constitute the main result of this paper. While F is not a strong function of lattice spacing 
around a ~ 0.5/ g^T {f3 ~ 8), it then turns over and falls roughly linearly in a at finer 
lattices. This indicates that the lattice spacing needed before the ASY scaling behavior sets 
in is around (3 = 12 {a = l/Sg'^T). 

If we insist on believing that the small a scaling behavior is of form ci + C2g^aT with 
C2 representing a correction to scaling, then a fit to the points with a < l/Sg'^T {(3 > 12) 
gives ci = 0.257 ± .044 and Ca = 3.29 ± .18 {^"^/u = 5.4/3). The "correction" to scaling only 
becomes subdominant below a = l/12.8g'^T {(3 > 51). The original motivation for believing 
in a finite small a limit for F is the belief that the UV lattice behavior is not important to 
the infrared dynamics. Such an enormous correction to scaling contradicts this belief. This 
makes it very difficult to reconcile our data with a finite small a limit for F. 

On the other hand, if ASY are correct, it makes more sense to plot our results as 
r/{a^aT^), as we do in Figure ^. The extrapolation to small a looks much better behaved 
here. If we fit F to the form F = a^aT^{ci + C2g'^aT) we get Ci = 74.5 ±3.4 and C2 = —73 ±12 
(^^/z/ = 1.6/3), which means the correction to scaling comes of order 1 at a = g^T {(3 = 4). 

In fact, as Bodeker has shown, we should not expect a finite intercept in this figure, there 
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should be a weak logarithmic divergence as a — > |§|. In the continuum its amplitude is ||T9 



log part of r = (10.7 ± .7)^ log (^] a'T^ . (13) 

The appearance of l/m'jj is changed somewhat on the lattice; m'j-) is reduced by a weighted 
average over k of the group velocity under lattice dispersion relations, which is about 0.68 
PJ. Combining this with the expression for m|, on the lattice, Eq. (0), gives 

log part of r = [7.8 ± .5] log ( -^-— J a^aT^ . (14) 

Strictly speaking such logarithmic behavior only pertains in the regime where log(l/5'^aT) ^ 
1, which is probably not satisfied at any conceivable lattice spacing. Nevertheless we perform 
the fit to see what happens. The fit in Fig. ^ which curves assumes F = log part of F + 
a^aT^{ci + C2g'^aT), which still has only two free fitting parameters. The best fit value is 
Ci = 54.6±3.4, C2 = — 39±12 {^'^/i' = 2.1/3). In this case the leading ci behavior dominates 
over the scaling correction out to a = lAg'^T (/5 = 3). Note that the data do not yet justify 
either believing or disbelieving in the log behavior. This is because the coefficient of the log 
is numerically small. 

4 Small volumes 



Ambj0rn and Krasnitz report another puzzling problem with the ASY picture |T3[. They 
analyzed the dependence of F on lattice spacing at a fixed, small volume, small enough that 
F is much smaller than its large volume limit. The idea was that, in such a small volume, iVcs 
fiuctuates about an integer value with occasional, abrupt changes from integer to integer, 
which can be identified even with the naive definition of A^cs- Therefore, A^cs can be tracked 
topologically without needing a true topological definition. Their results, replotted here in 
Figure |, are puzzling. Not only is there a strong lattice spacing dependence in evidence; 
it is too strong. The rate appears to fall off faster than F oc a, and certainly faster than 
the rate falls off in large volumes, where the corrections to scaling at finite a make the a 
dependence somewhat weaker than linear. 

To see why this result does not jive with the ASY picture, we will review again the basic 
ASY argument. Examining the propagator for the gauge field at momentum p gives an 
equation of motion for A{p) which, viewed on suitably long time scales, is approximately 

<fA{p) nml) dA{p) 



df^ ip dt 



-p'A{p) . (15) 



The damping term, with the single derivative, is more and more important as p gets smaller. 
We expect that, in a constrained volume, the field configurations responsible for changing 
A'cs are spatially smaller than in a large volume. Therefore they are composed of excitations 
of larger p, and should be more weakly damped. In particular it should take a larger value 
of m\) (smaller a) for the overdamped regime to apply, and the corrections to linear in a 
scaling should be larger. 
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Figure 5: Dependence of F on a in a small volume, L ~ A/g'^T. Left, the results Ambj0rn 



and Krasnitz reported in |jT3|. Right, results when more care is taken to keep the lattice 
volume fixed in physical units. The overly strong a dependence in the Ambj0rn-Krasnitz 
data is due to small changes in the physical volume, which F depends on very strongly. 



The problem lies in the data. Ambj0rn and Krasnitz worked in a lattice with iV = /3l, 
which at the unimproved level means L = A/g'^aT. However, they used the tree relation 
between the lattice spacing and the reciprocal temperature /3l- As can be seen from Eq. |12|, 
the true lattice spacing is larger than the tree relation indicates, and the error is worse as 
the lattice is made coarser. Therefore, their larger a lattices possessed more physical volume 
than their smaller a lattices. The effect is fairly small; the difference in linear dimension 
between their /?l = 10 and /3l = 20 lattices is only about 3%. However, as Fig. |] shows, F 
is a very strong function of volume in the regime where they were working. Around N = f] 
(L = A/g'^T), the data in the figure give roughly (i(logF)/(i(logL) ~ 10, so a 3% change in 
length could lead to a 30% change in F, which is significant. In fact this effect could be as 
large as or larger than the a dependence due to hard thermal loop dynamics. 

To fix this problem we recompute F in a fixed, small volume, but using the improved 
relation between the lattice spacing and the reciprocal temperature. We choose to use a 
volume equivalent (at the improved level) to the volume of Ambj0rn and Krasnitz' finest 
volume, /?L = 20 and A^ = 20, which, using the improved relation, gives L = A.llA/g'^T. 
The results are presented in Figure ^, and include a re-calculation of the finest lattice point, 
using a topological definition of Nq^ (rather than counting integer winding changes by eye 
from the old "unimproved" Nqs definition). The new data show a lattice spacing dependence 
similar to the large volume case, although to make a good comparison we would need data 
at smaller a (larger (3). The behavior of F in small volumes is in accord with what we expect 
if the ASY argument is correct. 
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5 Conclusion 

Chern-Simons number diffusion in pure Yang-Mills theory follows the Arnold-Son- Yaffe scal- 
ing behavior, F oc a; however it takes a fairly fine lattice to demonstrate this in a convincing 
way. The same behavior, with similar corrections to scaling, is observed in small volumes, 
but only after taking care to keep the physical volume constant beyond leading (tree) level 
while the lattice spacing is varied. 

The original goal of measuring F on the lattice was to determine the rate at which a 
baryon number excess, present before the electroweak phase transition, would be erased. We 
can do this with the results of this paper by using Arnold's study of the matching between 
the lattice and continuum theories p] . He shows that F on the lattice matches the continuum 
value when 

(.68 ± .20)m|i(latt) = m|, (continuum) , (16) 

where the error is all systematic, but the error estimate is considered generous P]. In the 
minimal standard model (MSM), at leading order mj^ = {ll/6)g'^T'^, which for g"^ = 0.42 
and using Eq. (0), means the MSM value is obtained on a lattice with a = .157 /g'^T, or 
f3 = 25.4. Our data at /5 = 24 give F = O.SSa^T^, with an error insignificant compared to 
the 30% error estimate in the lattice to continuum matching procedure, since F oc a is well 
satisfied in this regime, our estimate for the Standard Model value of F is 

F(MSM) = (0.82 ± .24)a^T^ or F = (45 ± 13) ( ^^ ) a^T\ (17) 



TTl 
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where the latter form shows the correct parametric dependence on the Debye mass. This 
result is in good agreement with results obtained when classical Yang-Mills theory is sup- 



plemented with particle degrees of freedom ||T0[ to induce the hard thermal loop effects. 



Therefore the diffusion constant for Ncs in classical Yang-Mills theory is consistent with 
both analytic expectations and numerical results obtained by explicit inclusion of hard ther- 
mal loops. 
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